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Abstract 

We explore the phase diagram of an 0(n) model on the honeycomb 
lattice with vacancies, using finite-size scaling and transfer-matrix 
methods. We make use of the loop representation of the 0(n) model, 
so that n is not restricted to positive integers. For low activities of the 
vacancies, we observe critical points of the known universality class. 
At high activities the transition becomes first order. For n = the 
model includes an exactly known theta point, used to describe a col- 
lapsing polymer in two dimensions. When we vary n from to 1, we 
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observe a tricritical point which interpolates between the universaUty 
classes of the theta point and the Ising tricritical point. 
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1 Introduction 



In two-dimensional models, exact solutions and Coulomb gas methods have 
led to the discovery of a large collection of universality classes, including some 
that depend on a continuously variable parameter (see e.g. Refs. such 
as the Baxter model, the critical and tricritical Potts models, and the critical 
0(n) model. However, such results are not known for the tricritical points of 
the 0(n) model with general n, in spite of the fact that there are no obvious 
reasons why such a family of tricritical points should be absent. Tricritical 
behavior is known only for discrete values of n: the so-called theta point at 
n = 10] and the Blume-Capel tricritical point 0, |§] at = 1, where the 
model intersects with the q = 2 Potts model, so that the results for Potts 
tricriticality apply. The branch of g-state Potts critical points applies to a 
limited range of q values. Its analytic continuation beyond the end point at 
g = 4 appears to cover the same range of q, and is then found to describe 
a branch of tricritical points In the case of the 0(n) model on the hon- 
eycomb lattice, the situation is similar in a sense: two branches of critical 
behavior are known [|^, which are the analytic continuation of each other, 
and connect at the upper boundary at n = 2. However, the analytic contin- 
uation of the critical branch does not describe a phase transition: instead, it 
describes the low-temperature 0(n) phase which is still critical in the sense 
that the correlations decay algebraically for general n. 

Several other branches of critical behavior are known for 0(n) models with 
additional interactions, see e.g. Ref. ^ for an analysis of the square-lattice 
0{n) model. One of these branches, called branch 3, marks the boundary 
between a range of ordinary 0(n) transitions and a first-order range, just as 
expected for tricritical points. However, this higher critical point of branch 
3 is accompanied by the freezing out of Ising-like degrees of freedom |§, |^ 
and does not match the known tricritical behavior at n = and 1. 

Another family of critical points, called branch 0, describes a dilute form 
of the 0{n) model. It includes a point at = that belongs to the univer- 
sality class of the theta point. In this sense branch is a generalization of 
the theta point to other values of n. Thus it could be a logical candidate for 
a class of tricritical points, but the n = 1 point of branch represents the F 
model [|lOl and not the tricritical Ising model. 

This somewhat puzzling situation raises some questions, such as what 
precisely is an 0{n) tricritical point, and what sort of connection exists be- 
tween the theta point and the Ising tricritical point? As for the first question. 
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we adhere to the accepted picture at n = and n = 1 that a gradual increase 
of the activity of the vacancies will eventually drive the 0(n) phase transition 
first order; the point where the first-order line meets the line of continuous 
transitions is considered a tricritical point. As for the second question, we 
may look for a regime of tricritical points that interpolates between a theta 
point and an Ising tricritical point. In the absence of exact results we use 
a numerical method, namely finite-size scaling of transfer-matrix results, in 
order to explore such tricritical points for general n. 

Thus we select a model for our study that includes vacancies explicitly. 
It includes the theta point defined by Duplantier and Saleur and the 
honeycomb 0(n) model without dilution as special cases. For n = 1, 2, ■ • • 
this model is equivalent with an 0{n) symmetric spin model |]TT|. The model 



and the numerical methods are described in Section 2. The results and their 
interpretation are presented in Section 3, and our conclusions and outlook 
are summarized in Section 4. 



2 The model and numerical procedures 

We define a dilute 0(n) loop model in which the loops run on the edges of 
the honeycomb lattice. Each edge covered by a loop has weight w, and each 
loop has a weight n. Dilution is represented by means of face variables: the 
elementary hexagons of the lattice are 'vacant' with weight y and 'occupied' 
with weight 1 — y. The loops are forbidden to touch any vacant faces. 

The occupied faces form a lattice C which is any subset of the dual (tri- 
angular) lattice. The loop configurations of this model are represented by 
graphs Q consisting of non-intersecting polygons, which avoid the edges of 
the vacant hexagons. Thus the partition function is 

z = j:j:y''^i^-yf'''^^''-^''' (1) 

c g\c 

where N^j is the number of empty faces specified by the vacancy configuration 
C The number A^^; of edges covered by Q is equal to the number of vertices 
visited by a polygon, and Ni the number of loops. The second summation is 
only over those graphs Q allowed by C 

For the construction of the transfer matrix we choose the usual geometry 
of a model wrapped on a cylinder, such that one of the lattice edge directions 
runs parallel to the axis of the cylinder. The finite size L of the model is 
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taken to be the number of hexagons spanning the cyhnder; the unit of length 
is thus a/3 times the nearest-neighbor distance. The transfer-matrix index 
represents the loop connectivity, i.e. the way in which the dangling bonds at 
the end of the cylinder are mutually connected, and the vacancy distribution 
on the last row of the lattice. The introduction of vacancies leads to an 
increase of the number of possible values of the transfer-matrix index for 
a given system size. Therefore the present calculations are restricted to 
somewhat smaller system sizes than e.g. those investigated in Ref. [^]. 

Our numerical analysis focuses in particular on the asymptotic long- 
distance behavior of the 0(n) spin-spin correlation function. In the language 
of the 0{n) loop model, this correlation function between two points sepa- 
rated by a distance r assumes the form g{r) = Z' /Z, where Z is defined by 
Eq. (|l|) and Z' by a similar sum on graphs Q' which include, apart from a 
number of closed loops, also a single segment which connects the two corre- 
lated points. We take these points far apart in the length direction of the 
cylinder. For the calculation of Z' we use the 'odd' transfer matrix built on 
connectivities that describe, in addition to a number of pairwise connected 
points, also a single point representing the position of the unpaired loop 
segment. The 'even' transfer matrix uses only connectivities representing 
pairwise coupled points. 

Let the largest eigenvalue of the even transfer matrix for a system of size 
L be and that of the corresponding odd transfer matrix A^^\ These 
eigenvalues determine the magnetic correlation length C,h{L) along a cylinder 
of size L: 

C^(L) = (2/v^)ln(Ai°VAi'^) (2) 
where the geometric factor 2/ a/3 is the ratio between the small diameter of an 
elementary face (a/3) and the length added to the cylinder by each transfer 
matrix multiplication (3/2). The present calculations of the relevant eigen- 
values of the even and odd transfer matrices are restricted to translationally 
invariant eigenstates. 

The magnetic correlation length could thus be calculated numerically for 
finite sizes L up to 10. Its dependence on the vacancy weight, the bond 
weight, and L allows the exploration of the phase diagram. Suppose that the 
parameters of the model are near a renormalization fixed point, at a distance 
t along some temperature-like field direction. The finite-size scaling behavior 
of the magnetic correlation length is 

C^(t,L)=L-ieV(tL^*) (3) 
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For sufficiently small t we may expand the scaling function 

4-^(t,L) = L-i(eV(0) + atL^* + ---) (4) 

where a is an unknown constant. For conformally invariant models, ^/^^(O) = 
2TTXfi where is the magnetic scaling dimension of the pertinent fixed point 
T^, O. This leads us to define the scaled gap Xh(t,L) = L/(27r^/i(t, L)) 



which scales as 

X,,{t,L)=Xh + aW^ + --- (5) 

where further contributions may be due to irrelevant scaling fields, or a second 
relevant temperature-like field if the system is close to a tricritical point. In 
actual calculations, we obtain the scaled gaps as a function of y and w, and 
their interpretation in terms of scaling fields involves one relevant field in the 
case of an ordinary critical point, and two in the case of a tricritical point. 
In order to get rid of the most relevant field we may e.g. fix y to some value 
and solve for w in 

Xh{y,w,L) = Xh{y,w,L + l) (6) 

The solutions, which are denoted w{y, L), are expected to approach the actual 
critical point as w{y, L) = w{y, oo) + cL^*'2~^* where Yt is the renormalization 
exponent of the most relevant temperature field and Yt 2 the second largest 
temperature-like exponent. Thus, by solving Eq. for a range of L values 
we obtain a series of estimates of the critical point w{y,L). In the case 
of an ordinary critical point, the irrelevance of other temperature-like fields 
(Yt^2 < 0) implies that the values Xh{y, w, L) taken at the solutions of Eq. (H) 
for large L converge to X^ as 

X^(t/,«;(y,L),L)-X,(xL^*'^ (7) 

By taking L) we have effectively removed the most relevant field. This 
formula still holds when y has a value such that y, w{y, oo) is in the vicinity 
of a tricritical point, but then the second temperature exponent is relevant: 
Yt^2 > 0. In that case, the values Xh{y, w, L) taken at the solutions of Eq. 
appear to diverge from the fixed-point value of Xh, i.e. the tricritical magnetic 
exponent. 

These considerations allow us to locate a tricritical point, if present. 
Convergent behavior of Xh{y,w{y, L), L) corresponds with ordinary criti- 
cal points, or with ffist-order transitions to which similar arguments apply. 
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Divergent behaviour, i.e. the Xh{y, w{y, L), L) are moving away from the tri- 
critical value of when L increases, reveals a tricritical point. Only when 
y is set precisely at its tricritical value, Xh{y,w{y,L),L) will converge. 

3 Numerical results 

A superficial exploration of the 0(0) model Eq. (|l]) by means of transfer- 
matrix calculations using small system sizes indicated that it has a phase 
diagram similar to that found in Ref. ||14[: the low-temperature and high- 
temperature phases are separated by a line of phase transitions which includes 
a range of ordinary 0(0) critical points, separated from a range of first- 
order transitions by an exactly known p, |^ theta point. Furthermore, this 
preliminary work indicated that the phase behavior does not change in a 
qualitative sense in the range < n < 1. Encouraged by this outcome we 
performed more extensive analyses for n = 0, |, and 1. We obtained solutions 
of Eq. (^ for fixed values y = 0, 0.1, ■ ■ ■, 0.8 using system sizes up to L = 10. 

The results for Xh{y,w{y, L), L) are shown in Figs. 1-3. For small y the 
data tend to converge to the known magnetic dimensions ^ for the 0(?7,) 
critical points, indicated by dashed lines. For large y values the data tend 
to converge to a much lower value of Xh which does not only reflect the 
expected behavior near a discontinuity fixed point [|15] but also the algebraic 
nature of the dense 0{n) phase. For intermediate y we observe the behavior 
formulated above for the vicinity of a tricritical point: the data indicate that 
there exists a special value of y for which the Xh{y,w{y,L),L) converge to 
a constant level which is different from the large- and small y asymptotes. 

For n = (see Fig. 1) the tricritical point is located ai y = |; the 
numerical data coincide with the exact tricritical value Xh = for all L. 
For n = I we do not know the exact location of the tricritical point and its 
associated magnetic exponent, but from the data in Fig. 2 we conclude that 
it is located near y = 0.54, w = 1.058, and that X^ ~ 0.036 (dotted line). 
For n = 1 we estimate the tricritical point at y = 0.58, w = 1.115. The 
finite-size data (see Fig. 3) are in a good agreement with the exactly known 
tricritical value Xh = 3/40 which is indicated by a dotted line. As expected, 
for increasing L the data are moving away from that line. 

The solutions w{y, L) seem to converge well with increasing L. The phase 
transitions thus found are shown as data points in the phase diagrams of Figs. 
4-6, which apply to n = 0, | and 1 respectively. The estimated errors are 
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much smaller than the symbol sizes. 

For y = our numerical procedures accurately reproduce the exact crit- 
ical points given in Ref. [|^: if = 0.54120 ■ ■ ■ for n = 0, w = 0.55687- ■ ■ for 
n = |, and w = 0.57735 ■ ■ ■ for n = 1. These exact critical points are shown 
as black circles. For n = (Fig. 4) the exactly known tricritical point (black 
square) coincides with the data point for y = ^. For n = ^ and n = 1 the 
estimated tricritical points are shown as black squares. The data points in 
each of Figs. 4-6 are connected by a line of phase transitions that divides 
into a continuous part (dashed line) and a first-order part (full curve). These 
curves do not represent any further data but serve only as a guide to the eye. 



4 Conclusion and outlook 

The numerical analysis of the 0{n) model with vacancies demonstrates that 
one can interpolate between the known tricritical points at = and n = 1. 
These two points are thus not isolated cases, but special cases of generic 0{n) 
tricriticality. In spite of the apparent difficulty to find a class of tricritical 
points analytically, the existence of such a class is obvious. Forthcoming 
numerical analysis will show whether this family of tricritical points covers 
the whole range —2 < n < 2 for which analytical results for the ordinary 
0(n) critical point are known. At present we see no obvious reason for a 
negative answer. 

The tricritical exponents remain to be determined; our present work con- 
tributes only Xh ~ 0.036 as a new result for the O(^) model. This determines 
the decay of the magnetic correlation function with distance r according to 
r"'' with 1] = 2Xh. As demonstrated e.g. in Ref. [§], several other exponents 
can be determined such as the temperature dimension Xt and magnetic expo- 
nents associated with different types of interfaces. A further characterization 
of the family of O (n) tricritical points will also require a determination of the 
conformal anomaly. While this is possible in principle |]16|, it remains to be 



seen whether the inaccuracies involved in the determination of the tricritical 
points will allow the emergence of a clear picture. The present calculations 
are still rather limited in magnitude; only a few hours of CPU time on a few 
Silicon Graphics workstations were used, which we found to be sufficient for 
the present purposes. Thus, the use of larger computer systems will open 
the possibility to explore significantly larger system sizes L, and thus help to 
conduct such more detailed investigations. 
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We conclude with a few remarks on the relation between 0(n) tricritical- 
ity and the percolation of the vacancies. For n = the vanishing loop weight 
enforces a perfect loop vacuum on the left hand side of the line of phase tran- 
sitions (Fig. 4). Thus, the site percolation on the triangular lattice apphes: 
the vacancy percolation line lies at ?/ = | and ends at the tricritical point. 
Once the density of the vacancies has increased above the percolation thresh- 
old, there is no way in which the loop configurations can reach criticality, and 
the transition becomes discontinuous. We thus expect that, also for n ^ 0, 
the vacancy percolation line connects to the tricritical point. However, for 
n > the loop vacuum is no longer perfect so that the percolation threshold 
moves to ?/ > |, in accordance with the estimated tricritical values of y. 
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Figure 1: Effective magnetic scaling dimension Xh of the 0(0) model ver- 
sus system size L, for various values of the vacancy weight y. The data 
points (not shown) are connected by lines. Counting from above, these lines 
correspond with y = 0, 0.1, 0.2, ■ ■ ■, 0.8. For small y the data behave in 
accordance with the expected magnetic dimension = 5/48 = 0.10416 • • • 
of the 0(0) critical point (dashed line). The data for y — 0.5 agree with the 
exactly known result for theta point of this model. 
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Figure 2: EfTective magnetic scaling dimension of the 0(i) model versus 
system size L, for various values of the vacancy weight y. The data points 
(not shown) are connected by lines. Counting from above, these lines corre- 
spond with y = 0, 0.1, 0.2, ■ ■ ■, 0.8. For small y the data behave in accordance 
with the expected magnetic dimension Xh = 0.11544- • • of the 0(|) criti- 
cal point (dashed hne). The tricritical magnetic exponent is estimated as 
Xh ^ 0.036 (dotted hne). 
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Figure 3: Effective magnetic scaling dimension of the 0(1) model ver- 
sus system size L, for various values of the vacancy weight y. The data 
points (not shown) are connected by lines. Counting from above, these lines 
correspond with y = 0, 0.1, 0.2, ■ ■ ■, 0.8. For small y the data behave in 
accordance with the expected magnetic dimension Xh = 1/8 of the Ising 
critical point (dashed line). The exactly known value of the tricritical Ising 
magnetic exponent X^ — 3/40 is indicated by a dotted line. 
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Figure 4: Phase diagram of the 0(0) model in the bond weight w versus 
vacancy weight y plane. The exactly known point at y = is shown as a 
black circle, that at y = ^ (the theta point) as a black square. The numerical 
data points are shown as open circles, connected by a curve. The dashed part 
indicates critical points in the 0(0) universality class, the full part indicates 
first-order transitions. 
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Figure 5: Phase diagram of the 0(|) model in the bond weight w versus 
vacancy weight y plane. The exactly known point at |/ = is shown as a 
black circle. The numerical data points are shown as open circles, connected 
by a curve. The dashed part indicates critical points in the 0(|) universality 
class, the full part indicates first-order transitions. The estimated location 
of the tricritical point is shown as a black square. 
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Figure 6: Phase diagram of the 0(1) model in the bond weight w versus 
vacancy weight y plane. The exactly known point at |/ = is shown as a 
black circle. The numerical data points are shown as open circles, connected 
by a curve. The dashed part indicates critical points in the 0(1) universality 
class, the full part indicates first-order transitions. The estimated location 
of the tricritical point is shown as a black square. 
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